Periodic Vortex Structures in Superfluid '^He-A 
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We discuss the general properties of periodic vortex arrangements in rotating superfluids. The 
different possible structures are classified according to the symmetry space-groups and the circula- 
tion number. We calculate numerically several types of vortex structures in superfluid '^He-A. The 
calculations are done in the Ginzburg-Landau region, but the method is applicable at all tempera- 
tures. A phase diagram of vortices is constructed in the plane formed by the magnetic field and the 
rotation velocity. The characteristics of the six equilibrium vortex solutions are discussed. One of 
these, the locked vortex 3, has not been considered in the literature before. The vortex sheet forms 
the equilibrium state of rotating ^He-A at rotation velocities exceeding 2.6 rad/s. The results are 
in qualitative agreement with experiments. 

PACS number: 67.57.Fg 



A superfluid cannot rotate homogeneously. Instead, 
quantized vortex lines are present in the equilibrium ro- 
tating state of superfluid ^He. In superfluid ■^Hc the ro- 
tating states are more-jdi verse. It has been discussed 
recently by Parts et. alH that four different types of vor- 
tices have been found experimentally in the superfluid A 
phase of "^He. In this paper we present theoretical studies 
concerning the vortices observed in '^He-A. 

Some of the theoretical results that we presented in 
Refs. El|2| were found to be incorrect in further calcula- 
tions. These errors are corrected here. As a consequence, 
the present phase diagram of vortices differs form the one 
in Refs. n^,0. In particular, there appears a new vortex 
structure, the locked vortex 3, but also the locations of 
other phase boundaries are changed. __ 

For Jntcoduction to superfluid •^HaJ'cl and its 
vorticegjij'El we refer to various review articles. Although 
not introductory, this paper intends to be a complete ex- 
position of what is needed for understanding the equilib- 
rium vortex structures in bulk superfluid "^He-A. 

We start in Section | with the formulation of the vortex 
problem, which is general enough for all superfluids and 
can be generalized also to superconductors. This gives 
a general classification of vortex states based on space- 
group symmetry and circulation number. The classifica- 
tion is continued in Section O using properties specific 
to ^He-A. The calculations of the vortex structures are 
based on the hydrostatic theory, which is discussed in 
Section |ll and the calculational method is described in 
Section tV, Detailed description of the different vortex 
types is given in Section |y|. The correspondence with 



experiments is discussed in Section VI 



I. THE GENERAL VORTEX PROBLEM 



Let us consider an uncharged fluid (in practice "^He 
or '^He) in a container rotating at angular velocity J7. 
We will neglect all complications arising from the finite 
size of the container. Although we will not discuss the 



detailed correspondence, the analysis in this section is 
also applicable to a charged fluid (superconductor) when 
Jl is replaced by the averaged magnetic field B. 

At the microscopic level, the fluid has the effective 
Hamiltonian HcS — Hq — il • J. Here Hq — X]i(-P?/2"i) + 
V is the Hamiltonian in a nonrotating system, which con- 
sists of a kinetic energy term and an interaction energy 
term V. The angular momentum J = X]J2(''* ^ P* ~ 
Pi X Ti) + Si] consists of an orbital and a spin part. We 
can write H^g in the form 
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where Vn^^ = $7 x r^ is the "normal fluid" velocity at the 
location of the particle i. The last term is the centrifugal 
energy. It causes the pressure to increase with increasing 
distance from the rotation axis. In principle, this term 
prohibits a strictly periodic vortex arrangement. How- 
ever, it is very small at experimentally relevant rotation 
velocities and container sizes, so that we can safely ne- 
glect it. We will neglect also the second-last term because 
it vanishes in '^He (S=0), and is very small in "^He, where 
it corresponds to a magnetic field of ~ 0.1 /zT at a typical 
H. = 1 rad/s. Because the rest of the paper is based on 
the reduced HcS, we write it again: 
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We will classify the rotating states according to their 
symmetries. For that purpose we first list all the sym- 
metries of the Hamiltonian (0). They are i) arbitrary 
translations, ii) arbitrary rotations around ft (ooz), iii) 
the combination (2^) of time inversion (') and rotation 
by the angle tt (2) around an axis perpendicular to J7, iv) 
the reflection (toz) in the plane perpendicular to H, v) 
the combination {m'^) of time inversion and reflection in 



a plane containing J7, and vi) all combinations of these 
operations. For each operation we have indicated a sym- 
bol in parenthesis. |-Ihey follow the international notation 
of crystallographyE3 with a prime added to denote time 
inversiontll. Throughout this paper we use a rectangular 
coordinate system xyz where the z axis is parallel to 17. 

It is not completely obvious that the translations per- 
pendicular to fl are symmetry operations of Hcs (y). 
This problem is equivalent to the case of electrons in uni- 
form magnetic fieldtj, and the corresponding problem for 
supcrfluid order parameter is discussed below. Another 
noteworthy feature is that both m^; and 2^ appear in 
combination with the time inversion. Otherwise these 
operations would not preserve the direction of the axial 
vector J7. 

Generally, the physical system either has all the sym- 
metries of the Hamiltonian, or alternatively, one or more 
of the symmetries are broken. An ordinary fluid would 
preserve all the symmetries of the rotating Hamiltonian 
(0). This is not the case for a superfluid. We will show 
below that at least part of the translation symmetry is 
broken in the superfluid state when fi 7^ 0. 

The fundamental property of superfluidity is that one 
quantum state becomes macroscopically occupied. This 
condensate is described by an order parameter A(r). The 
order parameter can be a scalar, as in *He, or a more 
complicated object. We associate a velocity field Vs to 
the particles in the condensate. There is no general ex- 
pression for the superfluid velocity Vg in terms of ^(r). 
Also, several different velocities Vg can exist, for exam- 
ple, one for spin up and another for spin down parti- 
cles. Irrespective of the precise definition, we only need 
to know how Vg changes in a gauge transformation. We 
require that the velocity associated with the order pa- 
rameter exp[i(/)(r)]A('')(r) is Vg = (S/M)V(/)(r) + vi°\ 

where Vg is the velocity corresponding to A^'^>(r). Here 
M is a mass that depends on the particular system. It 
equals the atomic mass for "^He {M — m^) and twice the 
atomic mass for ^He {M = 2TO3). 

It is now obvious that A{r) cannot be constant and 
also a phase factor exp[i(/)(r)] times a constant is not al- 
lowed. The reason is that the kinetic energy term in the 
Hamiltonian (0) would grow faster than linearly with the 
volume of the system because Vg would be constrained 
by f Vg • dr = and could not imitate Vn = Jl x r on a 
large scale. 

Our basic assumption is that the equilibrium structure 
of the rotating superfluid is periodic in space. It follows 
from above that the minimum period has to be finite at 
least in one direction, which is not parallel to ft. We 
do not make here any assumption whether the transla- 
tion symmetry is discrete or continuous in the two other 
linearly independent directions. 

Crystalline materials are classified according to their 
symmetry into 1651 magnetic space groups. The most 
effective way to label these is the international crystallo- 
graphic notationESIlil. We use the same notation to label 



the space groups of rotating superfluids. This is possi- 
ble because corresponding to every rotating state there 
exists at least one space group of a crystal. The reasons 
for this are that i) the symmetries of the rotating-fluid 
problem listed above are a subgroup of those possible for 
a crystal, and ii) no new symmetry groups appear even 
if one or two of the translation symmetries in a rotating 
superfluid were continuous. 

Not all the 1651 magnetic space groups are relevant 
for rotating fluids. Firstly, the time inversion is present 
in a rotating superfluid in a trivial way. In the generat- 
ing symmetry operations it appears in combinations with 
rrix and 2^, and only with them. Thus, by simply ignor- 
ing the time inversion, one can construct a one-to-one 
mapping from symmetry operations of a rotating fluid 
into sympaetry operations that do not contain the time 
inversiorO. Thus it is sufRcient to limit to the 230 crys- 
tallographic space groups, which do not include the time- 
inversion operation. Secondly, the number of relevant 
groups is further reduced because a rotation axis higher 
than 2 is allowed in the direction of Jl only. This implies 
that cubic groups are not acceptable. The remaining 194 
space groups each give rise to 1, 2, or 3 different symme- 
try groups of rotating fluids. This is because some of the 
crystal groups can be oriented in different ways relative 
to the J7 direction. 

What has been said above about symmetry does not 
directly apply to the order parameter A{r). The reason 
is that this complex quantity has phase (/>, which is not an 
observable quantity. Therefore, instead of being strictly 
periodic, A{r) is only quasiperiodic: 



A(r + afe) = exp[i0fc(r)]yl(r) 



(3) 



Here a^ are three linearly independent translation vec- 
tors (fc = 1, 2, and 3) and 0fe(r) are the corresponding 
phase shifts. Similar phase factors occur also in rota- 
tions, reflections, and time inversions. The quantity that 
has to be periodic in lattice translations is Vg — Vn. This 
gives a constraint for the phase shifts 0fe(r). Using the 
gauge invariance for Vg, one flnds V0fc(r) = Jl x a^. In 
order to simplify the formulas, we will repeatedly use the 
notation Ct ^(Af/fi,)0. The gradient of (/)fc is trivially 
integrated toliJ 



<?!)fc(r) = Cfe + X afc -r. 



(4) 



where Ck are constants of integration. An implicit re- 
quirement here is that Vg is defined on the path of inte- 
gration. We assume that the regions where Vg is unde- 
fined are at most one-dimensional. In this case it seems 
possible to choose the unit cell of the translation lattice 
so that Vg is well defined along all its edges. 

An important requirement is that the lattice- 
translation rule (|^) is consistent with a uniquely defined 
A{r). We express A{r -I- a^ -|- a^) as a function of A{r) 
using the translation rule twice. The result should be 
independent of the order in which the two translations 
by a^ and a^ are done. This gives the condition 
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where Nk are integers and Cijk the fuUy antisymmetric 
tensor. 

The lattice translation vectors a^ can be chosen in sev- 
eral different ways. Next we want to redefine the set 
{s-k} so that it is optimal for fm'ther analysis. The new 
as can always be chosen parallel to SI. Namely, setting 
aa = X;^ Nkaf"^, it follows from Eq. (|) that ag'x Q^O. 
Applying Eq. i^) to a new linearly independent set {afc}, 
we find that Ni = N2 = hut N3 ^ 0. The nonzero 
integer value of N3 implies that a continuous translation 
symmetry can exist only in the direction of CI. We can 
therefore additionally require that ai and a2 are primi- 
tive translation vectors, ie., they correspond to the min- 
imum (positive) value of il • ai x 3.2- 

We define the circulation number N as equal to N^ (H) 
corresponding to primitive ai and a2: 



ri • ai X a2 = ttN. 



(6) 



Similar to the symmetry groups, the different values of N 
can be used to classify the rotating states. N is called the 
circulation number because it is related to the circulation 
of the supcrfiuid velocity around a primitive cell 
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primitive cell 



where Vg — {M/h)vs. We note that the limitation to 
the boundary of the primitive cell in Eq. (R) is crucial in 
•^He-A, where the circulation is not generally quantized. 

Let us consider the case that there is a continuous 
translation symmetry along 17. This is an important case 
because all known vortex types belong to this category. 
However, very few general properties can be listed in ad- 
dition to those ones already mentioned above. The main 
simplification is that the primitive translation vectors ai 
and a2 can be chosen perpendicular to $7. These gener- 
ate a two dimensional Bravais lattice. Thus these rotat- 
ing states can be classified into five categories according 
to the symmetry of the 2-D latticeEJ: oblique, square, 
hexagonal, primitive rectangular and centered rectangu- 
lar. The number of possible space groups is considerably 
larger. In particular, the 17 two dimensional space groups 
listed in Ref. ^ are not sufficient for rotating states be- 
cause they lack the operation rriz . 

We illustrate the classification of vortices with known 
example cases. For a scalar order pararneter (^He) the 
Bravais lattice is hexagonal and N = lllj. It has the 
symmetry group P— — 7 — -, or shortly PQ/mm'm' . Gen- 
erally, the international svaabols consist of a letter fol- 
lowed by three symbol setaliJ. The letter shows the basis 
of the lattice, for example, P denotes a primitive and C 
a centered unit cell. The following three symbol sets de- 
scribe symmetries with respect to three different inequiv- 
alent axes, respectively. The first ^ or 6/to indicates 
that there is a 6-fold rotation symmetry and a reflection 



symmetry m, both with respect to the same axis. Here 
the 6 fold axis has to be parallel to Jl and thus the re- 
flection plane is perpendicular to O. The second set ^ 
describes a 2' symmetry and an m' symmetry with re- 
spect to an axis perpendicular to il. Finally, the third 
set ^ describes the same symmetries around the third 
inequivalent axis of the hexagonal lattice. 

The relative orientation of the space group and Jl is 
usually revealed by the primes because the primed axes 
are always perpendicular to fi. For some structures 
(C12'l, for example) pure symmetry considerations are 
insufficient to fix the direction of Jl. However, as proved 
below Eq. (P), the Jl axis always coincides with one di- 
rection of translation symmetry. This is not a conse- 
quence of symmetry but arises from the divergent rigid- 
ity of the voptex lattice against tilt deformation at long 
wave lengthalj. 

As another example, we consider vortices of ■^Hc-B. 
Two types of fully stable vortices are known. An iso- 
lated A-phase core vortex has symmetry class 00m' , and 
a double-core vortex 2m'm'\]\. Here the symmetry break- 
ing relative to — ^ arises from the core of each vortex. 
When these vortices form a lattice, the simplest possible 
structures have N = 1, and the space groups are Pdrn'm' 
and Cm'm'2, respectively. The lattice breaks the rotat- 
ing symmetry of the A-phase-core vortex to six-fold, and 
the two-fold rotation symmetry of the double-core vor- 
tex breaks the hexagonal lattice symmetry to centered 
rectangular. Both these effects are in practice very weak 
because the core sizes of the vortices in '^He-B are much 
smaller than the distance between vortices. 

The symmetry classification of vortices has previously 
been made only for point groups. This means that all the 
translations in the space group are ignored. Although 
this does not describe the whole lattice symmetry, there 
are several physical propertiosr&r which the point group 
gives a sufficient descriptionEtllS. We also comment on 
the notation. The symmetries 1, m' 2' often appear in 
dealing with vortices of ^He. Here 1 is the inversion and 
the sub-indexes in m'y and 2^ denote that they refer to the 
same axis. These operations give rise to five symmetry 
classes 1, T, m', 2', and ^. Here the first one contains 
only the unit element, the three middle ones have each 
one symmetry operation 1, m'y, or 2'y, respectively, and 
the last one has all the three (because two of them imply 
the third). In Ref. o, the same groups were labeled by 
letters o, u, v, w and uvw, respectively. Still another 
notation is due to Schonflies, and this was used to denote 
the same classes in Ref. |8[ Contrary to the international 
crystallographic symbols, these other notations do not 
allow a meaningful generalization to space groups. 

Let us study the meaning of the constants Ck in Eq. 
(^. We will show that Ci and C2 can be put to zero 
without losing generality. We consider an arbitrary order 
parameter fleld A(°)(r). We construct from it another 
held ^(r) by doing a translation by an arbitrary vector 
b as follows: 



^(r) = exp(iJi X b • r)yl(°)(r - b). 



(8) 



This field obviously has the same energy as the orig- 
inal one because the phase factor takes care that the 
counterflow velocity v = Vg — Vn is unchanged: v(r) — 
v(o)(r — b). By straightforward calculation one can verify 
that the coefficients Ci and C2 for the new field are re- 
lated to the old ones by C^ = Cj, ' +2flxh-ak. Choosing 
b appropriately, one can put Ci and C2 to zero. Thus 
the significance of Ci and C2 is that their values fix the 
position of the vortex solution relative to the rotation 
center. 

The coefficient C3 in Eq. (Q) is the phase shift in trans- 
lations parallel to fi. It is related to the superfluid ve- 
locity parallel to Jl. It often vanishes for symmetry rea- 
sons, but it can be nonzero for vortices of low symme- 
try. For example, consider a vortex with the symmetry 
group C12'l and a continuous translation symmetry in 
the z direction. The only point symmetry operation 2^ 
leaves the z component of the current unchanged. Thus, 
such a vortex generally has a nonzero net superfluid cur- 
rent in the z direction even though Vg^z = 0. Depending 
on the boundary conditions at z = ±00, this current 
may be compensated by a current arising from a nonzero 
Vs,z = Cz/az- 



II. SUPERFLUID ^He-A 



The continuous structures can be further classified by 
the numbers 
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These numbers are integers because d and 1 are peri- 
odic. They describe how many times the mapping from 
the primitive cell to the vectors d and 1 covers the unit 
spheres. 

The numbers N and vi are not independent. This fol- 
lows from the definition of the superfluid velocity 



Vs = y^m^Vr 



(12) 



(As above, we use v = (2m3/h)v, where 1113 is the mass 
of a ■^He atom.) This implies the Mermin-Ho relationEJ 



V X Vs = - ^ e^jkh^lj X VZft 



(13) 



ijk 



which together with Eqs. (|^) and ( [Til) gives 

N = 2i^i. (14) 



The previous section showed that the vortex structures 
in any superfluid can be classified according to the cir- 
culation number and the space group. In this section 
we continue the classification using properties specific to 
3He-A. 

The order parameter of bulk,superfluid ^He-A is a com- 
plex 3x3 matrix of the formHa 



A = Ad(m + in) . 



(9) 



Here d, m and n are unit vectors and m _L n. The 

amplitude A is a temperature and pressure dependent 
constant. It is conventional to define 1 = m x n, so that 
m, fi and 1 form an orthonormal set. 

As a first step, the vortices are classified to "continu- 
ous" and "singular" . The former alternative means that 
the bulk form (M) with constant A forms a good approx- 
imation to the order parameter everywhere in the primi- 
tive cell. The latter alternative means that this is not the 
case. This classification may not be precise in general, 
but there is no difficulty for the six vortex types to be 
considered here: only the "singular vortex" is singular, 
the other four are continuous. 

We note that only singular vortices exist for a scalar 
A because the amplitude of A has to vanish somewhere 
within the primitive cell of iV = 1. Continuous vortices 
are possible in ^He-A because nonzero circulation can be 
generated by appropriate m(r) and n(r) fields. 



III. HYDROSTATIC THEORY 

For a quantitative determination of the vortex struc- 
tures we use an energy functional F{A). In principle, it 
can be calculated from the effective Hamiltonian (||) as 
F{A) = -rin[Trexp(-i/cff/T)]. This is a functional of 
A{r) because the trace (Tr) is restricted to states hav- 
ing a given macroscopic A{r). Various approximations 
for F{A) are available: quasiclassical weak-coupling and 
weak-coupling-plus models, and phenomenological theo- 
ries such as the Ginzburg-Landau theory and the hydro- 
dynamic theory. 

The basic assumption of the hydrodynamic theory is 
that the deviation of the order parameter A from the bulk 
form M) is small. For this we have to require two condi- 
tions, (i) The magnetic field H should not be too large. 
In practice this condition excludes only a small region 
near the superfluid transition temperature. Tc, where the 
A phase is distorted towards the Ai phaseQ. (ii) The vec- 
tors d(r), m(r) and n(r) are sufficiently slowly varying 
functions of the location r. This implies that the hydro- 
dynamic approach can be used for continuous vortices, 
but it is insufficient for singular ones. 

Because of the slow variation, only terms up to the 
second order in the gradients of d, rh and n are needed 
in the energy functional. The functional can be written 
as 



^ = -^/^rfM/d + /h + /g). 



(15) 



Here the volume V of integration is assumed to consist 
of an (arbitrary) integral number of unit cells- The mag- 
netic dipole-dipole interaction /d is given bynd 



(16) 



(17) 



U = 2.9d|d X i|2. 
The magnetic anisotropy term is 

/h = igh(d-H)2, 
and the gradient energyn^ 

2/g = P±v2 + (p||-p^)(i.v)2 

+2Cv • V X i - 2Co(i • v)(i • V X i) 

+ /fb|ix(Vxi)|2+if5|(i.V)d|2 

+A-6^[(ixV).d,)f. (18) 

The gradient term includes also the kinetic energy, which 
is a function of the counterflow velocity v == Vg — Vn. 
It follows from the structure of the functional that 
Feq(H,fl) of the equilibrium state is a non-decreasing 
function of both H and il. The energy is normalized so 
that i^cq(H,0) = 0. 

It should be noted that Vg and 1 = rh x n are not com- 
pletely independent variables but are constrained by the 
Mermin-Ho relation (^3|). In order to avoid such com- 
plicated constraints, we use d, ih and n as the basic 
variables. With these variables the constraints are sim- 
pler: d, m and fi have to be unit vectors and rh _L n. 
The energy functional (Rq) can be expressed as a func- 
tion of d, rh and n. Only the gradient terms require some 
calculation, and we get 

2/g = (pil+2X7)^(m-An)' 

i 

+ {K, + Kr) ^[(m • Dm,)' + (n • Bfnf] 

i 
i 

+ {Kt + K^) ^[(n • T>m,f + (m • Dn,)'] 

i 

— 2Ki /^(n • T>rhi)m ■ Dn^ 

i 

+ (Kb + C - Co + Kr)[{B ■ m)2 + (D • n)^] 
+ 2{Kh ~Co + K^) Y^iCD ■ m)m,n • Bfn 

i 

— (D • ri)rhim ■ Dfii] 
+ (Co - C)^[{Dimk)Dkmi + {Dihk)Dkni] 

ik 



ik 

+ if5^(v,4)' 

ik 

+ {Ke - K5)J2ii^ ■ ^d,f + {h ■ Vd,f]. (19) 



Here we use gauge invariant derivatives Dm^ = Vm^ -|- 
Vafii and Dhi = 'Vfii ~ Vnmi. We have also used the 
notations Vn = {2m^/h)Q, x r, p = {h/2m3)'^p, C = 
(n/2m3)C, and Kt = p^_- p,, ^K,~ Kt + if b - 2Co. 

One can transform the gradient energy dl9) by par- 
tial integration. For example, / (i'^r(Vimfc)Vfcmi — 
J (firCV ■ rh)^ plus a surface term. For the present pur- 
poses such partial integrations can be done without pay- 
ing attention to the surface terms. The reasons are that 
i) the surface terms can affect the equilibrium configura- 
tion only near surfaces, if anywhere, and ii) although the 
local energy density is changed in the partial integration, 
all the energies of vortices are unchanged because of the 
periodic boundary conditions. 



The full gradient energy (19) is written down here 
in order to demonstrate that our calculational method 
is feasible whenever the hydrodynamic approximation is 
valid. In particular, the theory applies to all temper- 
atures T < Tc except a small region near Tc (due to 
the Ai phase) and another region around T = 0. How- 
ever, the present numerical calculations arc made in the 
Ginzburg-Landau (GL) regionllll. This means tempera- 
tures only near Tc (Tc — T <C Tc), but this range is still 
wider than the one that has to be excluded because of 
distortion towards the Ai phase. In this region the GL 
theory is more general than the hydrodynamic one. If one 
makes the hydrodynamic approximation in the GL the- 
ory, one arrives at the set of equations presented above, 
but with certain restrictions on the coefficients. They are 
P||/2 - Pi/(7 + l) = C/{j-ii-l) - Co/(7-l) =K,= 
Kt = Ki,/"f = K5/2 = Kg/{^ + 1). These conditions 
imply that all terms that are higher than second order in 
m or n disappear from the gradient energy (Il9). 

Collecting all simplifications and reducing units, we 
can write the energy terms as 

/d = i[(d.m)2 + (d.n)2] 

/h = ^(d-H)2 

4/g = ^[ViTOfe -I- (wn),nfc]' + ^[Vihk ~ Mimk]"^ 

ik ik 

+ (7 - 1)[(V • m + v„ • n)2 + (V • n - v„ • m)^ 
+ ^(m-Vd;)2-f^(n.Vrf,)'] 



+2 ^{Vidkf + 4(277 - l)m X n • J7, 



(20) 



where Vn — ftxr. We have used "dipole units" for length 
(Cd = (?i/2m3)i/p||/5d), field {H^ = ^gd/ffh), angular 



velocity (fid = fi/2m3^j), and energy density (gd). 

It should be noted that the functional (|2(]| ) is based 
on purely phenomenological considerations. These leave 
open only two dimensionless parameters 7 and 77. (In the 
notation of Ref. p2, 7 = K^/Kt and rj — Kc/Kt) They 
can be calculated in the quasiclassical theory. In the 
weak coupling approximation 7 = 3 and 77 = 1. A more 
complicated weak-coupling-plus model gives 7 w 3.1, but 
T] is unchangec£3. We have made a few tests that our 
results are not sensitive to the value of 7, so we will use 
the weak-coupling values in the following. 



IV. NUMERICAL CALCULATION 

For numerical computation the energy functional ( |20| ) 
was discretized on a lattice. We assume that there is no 
dependence on z, so that a two-dimensional square lattice 
(xj, yk) is sufficient. The first thing in the numerical pro- 
gram is to specify the magnetic field, the rotation veloc- 
ity f2, the circulation number N and the Bravais lattice. 
For rectangular lattices one needs to specify the ratio 
b/a, where a denotes the length of the shortest possible 
primitive vector, and b is the lattice constant in the per- 
pendicular direction in the rectangular cell. No oblique 
lattices were found. The area of the primitive cell was 
then determined from Eq. (m. For simplicity of bound- 
ary conditions, it was useful to choose the lattice constant 
of the calculational lattice commensurate with the prim- 
itive cell of the vortex lattice. The last preparatory step 
was to give an initial guess for the fields d, rh and n. 
Then the values d{xj,yk), m{xj,yk) and h(xj,yk) were 
changed iteratively. At each lattice point {xj,yk) in the 
primitive cell, a torque acting on the spin vector (d) and 
on the orbital vectors (rh and n) was calculated, aad the 
vectors were rotated proportionally to the torqueo. The 
proportionality coefficient was chosen experimentally to 
achieve fast convergence. The values of the fields out- 
side the primitive cell were determined from the periodic 
boundary conditions (||-Q). As discussed above, the co- 
efficients Ci and C2 are arbitrary but it was natural to 
choose them consistent with the initial guess in order 
to avoid unnecessary motion of the vortex. We assume 
C3 = 0. The iteration was continued until the energy 
converged, and the torques approached zero. 

For rectangular lattices, the energy _F needs to be min- 
imized also with respect to the ratio of the lattice con- 
stants u — b/a. This process is considerably simplified 
by noting that a calculation at a given value of u not 
only gives F{u), but also the first derivative dF{u)/du at 
constant area ab. It follows from the stationarity of the 
energy functional that 



dFiu) 
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2UV Jy 






Here Aij is the general order-parameter matrix in super- 
fluid '^He and c.c. the complex conjugate. Application to 
functional (EOh gives 



du uV Jv ^ 



(22) 



where 



4Ti = ^(V^TOfe + Vmfikf + ^(V^fife - Vmrhkf 



+ (7-1) 



{Viifii 



in,)(V 



+ (yini - Vnimi){V • n - Vn • rh) 

+ ^(mjVi(ife)rh • Vdfc + ^(n.,Vi(ifc)n • Vdk 

k k 

+2Y,i^^dkf■ (23) 

fe 

The iteration was often started with a rather small 
number of lattice points (~ 1000), and this was iterated 
to convergence. Then new lattice points were added in 
between the old ones, and the iteration was continued. 
The maximum final lattices contained around 100 000 
points. 

The procedure was repeated for several different initial 
guesses and values of magnetic field and rotation velocity. 
Although the finding of the minimum-energy structures is 
a well defined mathematical problem, physical intuition is 
needed in inventing the initial guesses. Particularly good 
guesses are the models used in previous investigations of 
vortices. In addition, we tried several variants of these. 
The initial guesses leading to the different structures are 
listed in the Appendix. 

V. RESULTS 

In all calculations T k, Tc and the field H was chosen 
parallel to the rotation axis z. The phase diagram of 
vortices in the plane formed by the magnetic field H and 
the rotation velocity Vt is shown in Fig. |l|. There are six 
equilibrium vortex types, which arc discussed separately 
below. In the names of the vortices, we follow Ref. [^. 



an A DyA,_j+ c.c. 
dDyAi, 



(21) 



The energy densities F (15) of the vortices are expressed 
in reduced units defined under Eq. (20). We choose the x 
axis parallel to the shortest primitive vector of the two- 
dimensional Bravais lattice. 

(i) The locked vortex 1 {INlXh&S-Lhc quantum num- 
bers TV = 4 and vi — Vd — 2L3'Ejt3. Its most distin- 
guishing feature is the square Bravais lattice. The space 
group is P^y-mJ^ '-"^ shortly PA/nb'm'. The symmetry 
operations of this as well as other space groups are listed 
in Ref. ^ 

Qualitatively the structure can be understood so that 
the first thing to minimize is the dipole energy (Rw. This 
gives a "locked" configuration where d = 1 everywhere. 
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FIG. 1. The phase diagram of vortices in the plane formed 
by the magnetic field H and the rotation velocity fJ || H. The 
calculation is done using the dimensionless units given below 
Eq. (bd). The solid lines denote phase boundaries between 
different continuous vortices (LVf , LV3, CUV, and VS). The 
LV2 should become the equilibrium locked vortex at very low 
fl < O.OOOSJld. The dashed lines denote the phase bound- 
aries of the singular vortex (SV) against the LV3 and the 
CUV. These are calculated using two values of the SV energy 
parameter: c = 2.6 (upper dashed line) and c = 3.1 (lower 
dashed line) . The real units for Q, and H are calculated using 
Qd ~ 120 rad/s and H^ = 2.0 mT, which are estimated for 
the pressure of 29 bar. 



If the field is zero, then the energy of this structure arises 
solely from the gradient terms (|l^) . These are minimized 
by a smooth distribution where the gradient of 1 has the 
same order of magnitude everywhere. Numerical calcu- 
lations show that this is achieved in a square lattice with 
A^ = 4, see Fig. |. 

The structure can be interpreted to consist of four el- 
ementary units. These units are called Mermin-Ho vor- 
tices because of a reseaabling structure first described in 
a cylindrical containeiO. The boundary of each unit can 
be defined hy Iz — 0. Two of the four units have a cir- 
cular distribution of 1 and /^ > (1 • il > 0) . The other 
two have a hyperbolic distribution and Iz < 0. It follows 
from Eq. (13) that each elementary vortex contributes a 
unity to A^ = 4. A finite axial field makes the cores of 
the Mermin-Ho vortices shrink because the field favors 
the orientation d _L z at the borders of the elementary 
vortices. 

At i7 = we find the energy F = 4:.72n°-^^'^ in the 
range il — 0.005 . . .0.04. In the case of perfect locking 
(d = 1), the energy would be strictly linear in fl. We find 
that this is only an approximation because with increas- 
ing n, the kinetic energy ( p^ is reduced at the expense 
of the dipole energy (|l^) . Our energy can be compared 
with Ref. |lj, where a variational calculation gives a 30% 
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FIG. 2. Locked vortex 1 (LVl) in zero field. In figures 0-n, 
the short lines denote the projection of the unit vector 1 to the 
X — y plane. For clarity of figure, the heads of the arrows are 
omitted. The component l^ is positive (1 ■ $7 > 0) in regions 
where (ix, ly) has circular appearance, and negative in regions 
where (l^, ly) is hyperbolic. The dots denote equivalent points 
in the periodic lattice. The LVl has the square Bravais lattice 
and the space group Pi/nb'm' . 



and a model based on isolated Mermin-Ho vortices a 6% 
larger value. 

A possible competitor to the LVl structure is a hexag- 
onal lattice with N = 2i>i = 2i'd = 2o. The ansatz form 
has 1 Ij J7 at the borders of the Wigner-Seitz cell. The 
space group is either P6m'm' for a radial and P62'2' for 
a circular distribution of 1. In agreement with previous 
authors, we find that this structure does not coripeapQud 
to the minimum energy at any values of fi and i;/ll3'OE2. 

(ii) The locked vortex 3 (LV3) has the same quantum 
numbers A^ = 4 and j// = j/tj = 2 as LVl. The main differ- 
ence is that the lattice structure is primitive rectangular 
rather than a square. The space group is Pjp-^-^, or 
shortly Pb'a'n. This structure has not been previously 
studied in the literature. 

The LV3 structure can be understood as a modification 
of the LVl structure. When the cores of the Mermin-Ho 
vortices shrink with increasing magnetic field, there re- 
mains a large bending of d « 1 outside of the cores. The 
gradient energy in this region can be reduced by rear- 
ranging the Mermin-Ho vortices. In LV3 the Mermin-Ho 
vortices form infinite chains, as visible in Fig. 0. The 
chains, or sheets consist of alternating circular and hy- 
perbolic units. Outside the sheets, the d « 1 fields are 
nearly constant and parallel (or antiparallel) to x. 

In contrast to all other vortex transitions, the change 
between LVl and LV3 seems to be of the second order. 
The lattice ratio h/ a grows continuously from unity at 
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FIG. 3. Locked vortex 3 (LV3) for H = 0.3ffd and 
Q, = O.OOefid. The notation is the same as in Fig. ||. The LV2 
has the primitive rectangular Bravais lattice and the space 
group Ph' an. The spacing of the plotted 1 vectors is l.OS^d. 



the transition. At constant H, the b/ a ratio grows with 
decreasing Q, and reaches 3.4 at our lowest Q. = 0.001 
■At H = 0.6. We fit F = 0(3.62 + Q.i7n~^-^^) a.t H = 
0.3 in the range Q = 0.001 . . .0.008. Here all the three 
parameters are free in the fit. 

It was demonstrated by Fmita and Ohmi that the LVl 
is unstable to a deformationE^. They found a structure 
where four Mermin-Ho vortices form a unit that is sep- 
arated by some distance from the other units. We find 
that this structure has higher enexgy than the LV3. 

(iii) The locked vortex 2 (LV2)E3 is defined by iV = 2 
and 1/1 = 1/^ = 1- The space group is C12'l, or shortly 
C2', where C denotes the centering of the rectangular 
lattice. The LV2 represent an alternative deformation of 
the LVl when the cores of the elementary vortices shrink 
in increasing magnetic field. Here the Mermin-Ho vor- 
tices form pairs. Each pair consists of one circular and 
one hyperbolic unit, see Fig. m Outside such a pair, 
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FIG. 4. Locked vortex 2 (LV2) for H = 0.4ffd and 
n = 0.002r2d. The notation is the same as in Fig. H. The 
LV2 has the centered rectangular Bravais lattice and the space 
group a! . The spacing of the plotted 1 vectors is 1.97^d. 



the d « 1 fields are nearly constant and parallel to x. 
This special direction breaks the hexagonal lattice sym- 
metry, which otherwise could apply to such well sepa- 
rated vortices. Therefore, we expect the lattice struc- 
ture is centered rectangular. We find that the ratio of 
the two lattice constants hja w 1.8. This is less than 
hja = ^/6, which is obtained from the hexagonal lattice 
(having b/a = y/S) by scaling x and y according to the 
anisotropy of the superfluid density (p± = 2pu). 

The LV2 has only one point symmetry 2' which means 
a rotation by tt around the y axis combined with time in- 
version. In previous '^He literature the 2' symmetry was 
denoted by wB. The symmetry implies that dy and ly 
change signs when a; ^ — a; in Fig. H, and other com- 
ponents of d and 1 remain unchanged. The symmetry 
transformation of ih and n depends on the specific choice 
of the phase factors. 

The LV2 is doubly degenerate. The degenerate forms 
are obtained from each other by exchanging the posi- 
tions of the circular and hyperbolic Mermin-Ho vortices. 
Formally this can be done by the operation m'. In the 
calculation we have only studied the simplest case where 
all LV2 vortices have the same orientation. If both de- 
generate forms are present simultaneously, it would lead 
to a larger primitive cell or, in an extreme case, absence 
of periodicity. 

We find the energy F = fl[~l.U + 1.371n(l/f7 + 
1/0.00412)] at H = 0.3 in the range i} = 0.001 . . .0.008. 
In this fit we kept the constant 1.37 multiplying the log- 
arithm fixed. This constant arises from the flow far from 
the vortex cores, and it should become exact in the limit 
51 ^ 0. The numerical value 1.37 is calculated in Ref. 
|3l| . The constant -1.14 can be interpreted as the core en- 
ergy of the vortex line. The constant 0.00412 describes 
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FIG. 5. Continuous unlocked vortex (CUV) for 

H — O.Siifd and Q — 0.02f2d- The notation is the same as 
in Fig. pi Similar to LV2, the CUV has the centered rectan- 
gular Bravais lattice and the space group C2'. The spacing 
of plotted 1 vectors is 0.59^d. 



the angular velocity where the cores of the vortices start 
to overlap. In the range of our calculation the LV2 has 
higher energy than LV3. However, it is expected that 
the LV2 becomes absolutely stable in the limit 17 ^ 0. 
Extrapolating the expressions of the energies we can es- 
timate that this takes place at 17 < 0.0005. 



(iv) The continuous-uiilp 
lyi = 1, and i^a = OBm. 



Ltex (CUV) has iV = 2, 
3. The space group is 



C2'. The CUV is similar to the LV2 with respect to N, 
vi, and the space group, see Fig. |[ The crucial difference 
compared to the LV2 is that Vd vanishes. Vortices where 
d and 1 differ essentially from each other [yd ^ v\) are 
called "unlocked" . 

The principle in all unlocked vortices is to minimize the 
field energy (|7|) in the first place, and therefore they are 
more economical in large fields H > H^ than the locked 
structures. In the CUV, d is approximately constant and 
parallel to x everywhere. The region where d and 1 differ 
from each other is called a "soft core". Its size ~ ^d is 
determined by the balance of the kinetic energy of v ( |l^ ) 
outside of the soft core and the dipole energy ( p^ inside. 
(Note that the total gradient energy of 1 in the core is 
approximately independent of the size.) 

The first suggestion for CUV had the m' symmetry 
Similar to several previous calculations, we find that the 
2' symmetric, fcajjHpias a lower free energy at all values 
of Vl and iK3'L3E3c3. The m' symmetric form seems to 
correspond to a saddle point of the free energy. 

Similar to the LV2, the CUV probably has the cen- 
tered rectangular lattice, and we find that b/a ~ 2. 
The energy of the CUV depends surprisingly much on 



the magnetic field. At iJ = we fit F = ^[-0.728 -I- 
1.371n(l/17 -K 1/0.219)] but at H = 0.6 we find F = 
O[-0.591 + 1.37 ln(l/rj + 1/0.0894)] . The expressions are 
most accurate in the range ft = 0.008 .. .0.03. Similar 
to the LV2, the value 1.37 is kept fixed in fitting the two 
other parameters. The interpretation of the other param- 
eters is the same as for LV2. The curving of the phase 
boundary between CUV and VS (Fig. |^) arises from the 
field dependence of CUV, as VS seems to be rather in- 
sensitive to H. 

According to our earlier calculation, the transition be- 
tween the CUVj-and the LV2 takes place aX H — 0.4 in 
the limit 17 -^ 0c3. This agrees with the present calcula- 
tion. We find the triple point between the LVl, the LV3, 
and the CUV at H = 0.29 and 17 = 0.014. 

(v) lChe|_yortex sheet (VS) has TV = 4, t/; = 2, and 
Vd = GolIo. The likely space group is Ph'a'n, similar to 
the LV3. Here N, I'l, and Vd are all the same as for the 
CUV except multiplication by 2. The crucial difference 
between the VS and the CUV becomes evident when their 
primitive cells are stacked one after another: in the CUV 
the soft cores form a two-dimensional lattice of lines while 
in the VS they form a series of equidistant planes parallel 
to X, see Fig. g. The d vector is approximately constant 
and parallel to x everywhere. 

The close similarity of the CUV and the VS is illu- 
minated if one thinks bending a VS and closing it to a 
cylinder. The CUV represents the smallest among such 
cylinders because it contains just one periodic unit of 
a vortex sheet. This relationship is also evident in the 
ansatz forms given in the Appendix. 

We have discussed above how the LVl is transformed 
to the LV2 by pairing the Mermin-Ho units. Then the 
CUV evolved from the LV2 when the dipole locking was 
removed. After that the VS was developed by opening 
the cylindrical structure of the CUV. We can now return 
to the starting point by noting that forcing dipole locking 
in the VS gives the structure of LVl and LV3. The close 
similarity of the VS and the LV3 is evident from the 
symmetry groups and from Figs. and S. 

The 1 vector in the VS is approximately parallel to x 
outside of the soft cores, but it has opposite directions 
on the two sides of the sheet. This implies that the prim- 
itive cell must contain (at least) two neighboring sheets. 
This is the reason for the double size of the primitive cell 
compared to the simplest lattice structure of the CUV. 
The surface tension of the sheet makes the primitive cell 
rather short in the x direction. We note that the present 
definition of b according to the primitive rectangular Bra- 
vais lattice is twice as large as in Refs. |[^ and^, where 
b denotes the distance between neighboring sheets. Be- 
cause the sheets are far from each other, the sliding of 
the sheets relative to each other probably leads to a neg- 
ligible change in the energy. 

The properties of the VS depend very weakly on 
the magnetic field in the studied region H < 0.6. 
The results of our calculations can be represented by 
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FIG. 6. Vortex sheet ( VS) for i/ = 0.2//d and fi = O-OSf^d . 
The notation is the same as in Fig. eI The VS has the prim- 
itive rectangular Bravais lattice and the space group Ph' an. 
The primitive cell contains two vortex sheets parallel to x. 
The arrows denote the opposite directions of I on different 
sides of the sheets. The spacing of plotted I vectors is 0.39^d. 



hja = 1.26r2-"-255 and F = X.ihVi^-^^^ for Q. = 
0.016 .. .0.07. Both these quantities are slightly larger 
than obtained from the twist-section modelcl, which gives 
hja = i(18/r2)i/3 and F = i(18r22)i/3_ The uniform 



windinaanodel is a slightly more complicated variational 
ansatzQEj. It gives an upper bound for the energy that 
is 10% higher at f7 = 0.016 and 3% higher at fi = 0.07. 
In the limit iJ ^ 1 we find the transition between CUV 
and VS at f7 == 0.022. 

A new feature in the phase diagram is that both the 
CUV and the VS are stable also in zero field. Although 
it has not been stated explicitly, the crossing of the en- 
ergies of the LVl and the CUV at if = appears also in 
previous literature. Comparison of the energies given in 
Refs. 14 and 40, for example, gives it at fi = 0.023. It 
was calculated in Ref. ^ that there is a transition from 
the LVl to a singular z vortex when H = 0.11. This pre- 
diction has to be revised because both the CUV and the 
VS have a much lower energy at this rotation velocity. If 
there is a transition to the z vortex, it takes place at a 
much higher H than expected in Ref. ^. 

The energy difference of the CUV and the VS rela- 
tive to the LVl at if = arises from competition of the 
dipole energy and the gradient energy of d. The LVl 
is stable at a low i7 because the dipole energy is min- 
imized in the first place, and the gradient energy asso- 
ciated with d is not important. With increasing $1 the 
gradient energy becomes larger. At the transition point, 
it becomes more economic to arrange d approximately 
constant although it means increased dipole energy in 
the soft core of the CUV or the VS. Based on purely di- 
mensional considerations, this transition w£is-expected at 
r2 ~ 1 (rj ~ fid ~ 120 rad/s in real units)OJ. However, 
the present calculation gives the transition between the 
CUV and the LVl at fi = 0.019, which is almost two 
orders of magnitude smaller than the naive expectation. 

(vi) The singular vortex (SVl-liaS|_^ = 1 and the space 
group Clm'l, or shortly Cm'HLjO. No other vortex 
considered here has TV = 1 because it is not possible 
for the continuous vortices as a result of Eq. (|lj) . The fl 
dependence of the energy can be written F = n[c{H, T)~ 
0.70 In fi] at small fl. Here the factop-p.70 arises from 
the flow field far from the vortex lineE^I. This factor is 
approximately one half of the value for a TV = 2 vortex 
line. Therefore, the singular vortex is favored over the 
CUV at a low r2. 

Based on topological arguments alone, the structure of 
the singular vortex could be very simple. For example, 
rh + in = exp(i(/))(y + iz) and d = 1 everywhere except 
at the singular "hard core", where these quantities are 
not defined. Here (f> is the azimuthal angle. However, 
energetics prefers a more complicated jStructure that has 
a soft core in addition to the hard coreEHl. This is because 
a structure with radial 1 and constant rh = z has lower 
energy than the simple vortex. Outside of the soft core, 
d and 1 are both nearly parallel to x. Similar to the LV2 
and the CUV, this gives rise to a centered rectangular 
lattice. 

The hydrodynamic approximation used for the calcu- 
lation of the continuous vortices is insufficient in the hard 
core of the SV. Therefore we have not calculated the 
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FIG. 7. Sketch of the core of the singular vortex (SV). 
The arrows denote 1 and d. The soft core, where d and 1 
differ considerably, appears as shaded. 



m' which in previous '^He literature was called v. In 
terms of the vectors this means that dy and ly change 
signs when y —J — y, and other components of d and 1 
remain unchanged. The SV is doubly degenerate. The 
two forms are obtained from each other by the symmetry 
operation 2y. 

All the discussion above was for the case where the 
field H is parallel to il. As far as we knpsM— other di- 
rections are considered only for the CUVtJcJ and the 
VSe3. Generally, it can be expected that the effect of 
field direction is not large in unlocked vortices, where d 
is approximately constant, and thus /h lR7i) can equally 
be minimized for arbitrary direction of H. 



VI. COMPARISON WITH EXPERIMENTS 



function c{H, T) in the energy F. Contrary to the con- 
tinuous vortices, there is also a temperature dependence 
c{H,T) « c{H) + 0.71n(l - T/T^) because the size of 
the hard core depends on T. In order to get an idea of 
the complete phase diagram, the phase boundary of the 
SV is included in Fig. m (dashed lines) by two arbitrarily 
chosen constant values of c{H, T) . 

The core structure of the SV can in principle be cal- 
culated using the Ginzburg-Landau theory, but there are 
two difficulties. Firstly, this theory introduces additional 
parameters (such as the coefficients of the five energy 
terms that are of the fourth order in the order param- 
eter matrix) whose values are not well known. So the 
accuracy of the results would be less than in the hydro- 
dynamic theory of continuous vortices. Secondly, the nu- 
merical calculation is difficult because the length scales 
associated with the soft (fd ~ 10 /^m) and hard cores 
(« 10 nm) are very different. 

We have made numerical simulations with Ginzburg- 
Landau theory where the difference in the soft and hard- 
core scales is arbitrarily reduced. We cannot expect any 
quantitative results from such a calculation, but we be- 
lieve that the following qualitative results are valid inde- 
pendently of our approximation. The structure around 
the hard core is 



± 1 = y sin (/> -I- cos 0(x cos ry -t- z sin rf) , 



(24) 



where 7y is a constant pigle, see Fig. M for illustration. 
The original suggestions has rj = -k 12 whereas the min- 
imum energy of the structure (p3) prefers 77 = CO. The 
true structure is likely to fall between these limits be- 
cause ?7 = would imply a large gradient energy (18) in 



matching the hard core (24) with the constant 1 outside 
of the soft core. 

There is no circulation around the hard core. Thus all 
the circulation arises from the soft core, which is quali- 
tatively described by the numbers vi = 1/2 and Vd = Q. 
The centers of the hard and soft cores are displaced from 
each other in the x direction. The vortex has symmetry 



In order to compare the calculated phase diagram with 
measurements, one needs to know how to prepare the 
equilibrium state in experiments. This is not simple be- 
cause the energy barriers separating the different vortex 
types are genecally so large that it is difficult to induce 
any transitional. So the decay of a metastable vortex 
type to the equilibrium type may be so slow that it can- 
not be observed. Also, if the rotation is started in the 
superfluid state, the vortex type that nucleates is gener- 
ally not the equilibrium one. For example, only continu- 
ous vortices are nucleated if the rotation is started in the 
superfluid state; no singular vortex has been observed by 
this method. 

The only exception to the above seems to be the re- 
gion very near the superfluid transition temperature T^. 
There the energy barriers separating the different vortex 
types are smallest. A practical way to perform the exper- 
iment is to cool slowly from the normal state [T > Tc) to 
the superfluid state at constant fJ and H . This procedure 
is expected to yield a state near the equilibrium one. It 
is important to remember, however, that the details of 
the transition in the presence of a thermal gjcadient and 
a magnetic field may be rather complicatedtl. Another 
limitation of the experiments is that they do not resolve 
the difference between the three types of the LV. 

Also needed for the comparison are the values of f^d 
and H^. We estimate fid ~ 120 rad/s, and H^ « 2.0 mT 
at 29 bar pressure. These are based on a weak-coupling 
analysis corrected by the enhancement of the energy gap 
according to Ref. E3. In addition, we have used the mea- 
suredjsiiift of the transverse NMR resonance frequency in 
the Ae3 or the B phasecJ, both of which give essentially 
the same result. The gap enhancement, which in Ref. ^ 
is given for the B phase as a function of the speciflc heat 
jump ACb, is applied to the A phase by replacing ACb 
byfACA. 

The comparison of the experimental and theoretical 
phase diagrams is shown in Fig. 3 of Ref H. The quali- 
tative agreement is good. A slight difference is that the 
vortex sheet is not observed in the experiments at the 
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maximal angular velocity 3 rad/s although according to 
the present calculation it should show up above 2.6 rad/s, 
assuming i^d = 120 rad/s. Possible explanations are that 
rid is larger than we estimated or the cooling through Tc 
does not accurately produce the equilibrium state. 

The six vortex types discussed above seem to be able 
to explain all the experiments that have been made in 
rotating bulk "^He-A. We comment, here on one contro- 
versial experiment. Torizuka et al£3 observe a transition 
in the rotating state at fi = 3 rad/s when the rotation 
velocity was varied at constant H = 0. The original in- 
terpretation in the same reference postulated a layer of 
vortices on the container wall. We consider this interpre- 
tation unlikely because such a layer is probably unstable, 
as pointed out in Ref. |4^. The present calculation now 
gives the possibility that the transition could be from 
the LVl to either the CUV or the VS. Unfortunately, the 
collected experimental data doea_nDt seem sufficient to 
identify the structure at large fcZrH. 



VII. CONCLUSION 

We have presented numerical calculations of the vor- 
tex structures in the Ginzburg-Landau region, and con- 
structed a phase diagram in the H — Q plane. There is in 
principle no difficulty in extending these calculations to 
lower temperatures. Although the phase diagram is not 
accessible experimentally at low temperatures, the cal- 
culation of the vortex structure would form the basis for 
a calculation of the NMR frequency shifts, which have 
been measured accurately. 

Our search of vortex types was based on previous sug- 
gestions. There may well be structures which could not 
have evolved from the initial guesses we have used. In 
particular, only the simplest periodic structures were 
tested. From the experimental point of view, it seems 
that there is at present no need to introduce new types of 
vortices that are stable in bulk ^He-A. That may change, 
however, when new regions are studied and more accu- 
rate measurements are done. In particular, low temper- 
atures, high rotation velocities, high magnetic fields, the 
neighborhood of the Ai phase, and restricted geometries 
could be studied. The studies could also be extended to 
metastable structures. For example, the LV was identi- 
fied from its metastable modification in high field, which 
more appropriately should be classified as a new type of 
vortexH. 
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APPENDIX 

We give approximate expressions for the order param- 
eter in different vortex structures. These can be used as 
initial guesses to produce the stable vortices discussed in 
Section 0. 

All structures can be simply represented using Euler 
angles (a, (3, 7) but choosing x as the polar direction: 

1 = xcos/3-1- sin/?(ycosa -I- zsina) (25) 

m -I- in = [—X sin /3 + cos (3{y cos a + z sin a) 

+i(— y sina + zcoso;)] exp(— 17). (26) 

In these coordinates the superfluid velocity ( p^ ) is 

Vs = -V7-cos/3Va. (27) 

Depending on the vortex type, we present the Euler an- 
gles as functions of either rectangular {x, y, z) or cylin- 
drical coordinates (r, 0, z). Unless explicitly stated oth- 
erwise, the ansatz forms have the same symmetries as the 
vortices that they represent. 

The VS has d = x and a — —7 — 7r/2 — 2'kx sgn(y)/a, 
where sgn(j/) denotes the sign of y. (3 is a, monotonic 
function of y so that (3{—b/2) = — tt, /3(0) = 0, and 
/3(fe/2) == TT. Especially at low velocities this function 
is strongly nonlinear so that all the change of (3 takes 
place in narrow regions (thickness ~ 1) at the two vortex 
sheets, which are located at y — zLb/A. This form of the 
order parameter is for a gauge where Vn = —2Clyx. The 
transformation to the more usual gauge Vn = J7 x r is ob- 
tained by including an extra factor exp(iQxy) multiplying 
the right hand side of Eq. (Gq) . The lattice constants are 
constrained by ab = Air/n ffl. 

An approximation to LVl and LV3 is the same as for 
VS except that d = 1. The best guess for LVl corre- 
sponds to a = 6 and nearly linear P{y). In spite of these 
choices, this ansatz has less symmetry than the converged 
solution for LVl. A symmetric but more complicated 
ansatz was suggested in Ref. 14 , A third alternative is to 
glue together ansatzes of four Mermin-Ho vortices. 

The CUV has d = x and a — —7 = (f>. /3 is a mono- 
tonic function of the radius r, which has /3(0) = and 
/3(r) « TT for r > tq. Here the radius tq is of the order of 
a few units of length (^d)- 

The LV2 is similar to CUV except that d = L rg is of 
the order of H^^ in dimensionless units. 

The SV can be generated by d = x, 7 = 0, and 
a = n/2 + (j). /? is a monotonic function of r, which 
has /3(0) = 7r/2 and /3(r) « tt for r > tq. Here tq is ap- 
proximately unity. This ansatz form has more symmetry 
than is present in a converged solution. 
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